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Abstract 

An  activated-carbon  electric-double-layer  capacitor  is  characterized  and  an  adaptive  algorithm  is  derived  and  implemented  that  is  com¬ 
mensurate  with  potential  traction  applications  of  these  energy  storage  devices.  The  electrochemical  characterization  relies  on  a  simplified 
equivalent  circuit  interpretation  extracted  from  a  more-complete  mathematical  representation  of  the  capacitor  system.  In  addition  to  the  high 
power  capability  and  potentially  low  costs  of  this  class  of  capacitors,  we  clarify  the  substantial  invariance  of  the  device  performance  with 
respect  to  temperature,  a  distinct  advantage  over  battery  systems.  The  robustness  of  the  equivalent  circuit  in  terms  of  capturing  the  salient 
features  of  the  experimental  data  over  the  temperature  and  voltage  range  of  interest  enables  the  formulation  of  a  model-reference  adaptive 
algorithm.  The  algorithm  developed  and  implemented  in  this  work  is  fully  recursive  in  that  the  only  variables  required  for  on-line  regression 
are  those  of  the  previous  time  step  and  the  current  time  step.  Successful  comparisons  of  the  algorithm’s  adapted  state  of  charge  and  power 
predictions  provide  an  initial  validation  of  the  algorithm. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Although  activated-carbon  electric-double-layer  capaci¬ 
tors  yield  high  power  over  a  broad  temperature  range  and 
consist  of  potentially  low-cost  materials,  they  have  not  been 
characterized  sufficiently  from  the  viewpoint  of  automotive 
traction  applications,  and  the  control  of  these  systems  has 
not  been  studied  in  the  context  of  formal  control  algorithms 
(i.e.,  optimal  estimators  and  adaptive  filters).  The  motivation 
behind  determining  a  battery’s  (i)  state  of  charge  (SOC),  (ii) 
state  of  health  (SOH),  and  (iii)  power  capability  (PC)  for  both 
discharge  and  charge  has  been  described  recently  [1].  The 
same  motivation  applies  to  electric-double-layer-based  ca¬ 
pacitors.  The  SOC  corresponds  to  the  stored  charge  available 
to  do  work  relative  to  that  available  after  the  electrochemical 
system  has  been  fully  charged  and  can  be  viewed  as  a  thermo¬ 
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dynamic  quantity.  SOH  refers  to  how  well  the  system  is  func¬ 
tioning  relative  to  its  nominal  (rated)  and  end  (failed)  states. 
The  same  three  characteristics  (SOC,  SOH,  and  PC)  are  re¬ 
quired  for  the  rational  integration  of  electrochemical-based 
energy  storage  devices.  A  general  survey  of  electrochemical 
capacitors  can  be  found  in  Conway’s  text  [2].  In  this  work  we 
are  specifically  interested  in  capacitors  based  on  high  surface 
area,  activated-carbon  electrodes.  Helpful  reviews  concern¬ 
ing  the  utility  and  characterization  of  these  systems  can  be 
found  in  [3-5].  In  the  2002  paper  by  Chu  and  Braatz  [4], 
the  SAFT  carbon-based  capacitor  was  shown  to  be  the  best 
performing  of  the  non-aqueous-solvent  systems  investigated 
by  the  authors,  as  shown  clearly  in  the  Ragone  plot  (specific 
energy  plotted  against  specific  power)  in  Fig.  3  of  their  paper. 
The  SAFT  capacitor  is  further  characterized  in  this  commu¬ 
nication.  Our  goal  is  to  understand  sufficiently  the  under¬ 
lying  physical  chemistry  of  these  devices  so  as  to  allow  for 
vehicle  integration  and  control.  Related  to  automotive  propul- 
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sion  systems,  most  texts  covering  control  theory  will  treat 
motor  control  and  feedback  systems  associated  with  vehicle 
speed  control,  little  attention  has  been  given  to  the  adaptive 
control  of  battery  systems  [1,6,7])  or  electric-double-layer 
capacitors  [8]. 

We  have  chosen  to  investigate  this  class  of  capacitors 
because  they  show  potential  for  providing  a  much  higher 
specific  power  (kWkg-1)  and  power  density  (kWl-1)  than 
advanced  batteries  at  a  much  lower  cost  ($  per  kW).  The 
material  costs  should  be  much  lower  than  that  of  advanced 
batteries,  as  the  primary  costs  corresponds  to  activated  car¬ 
bon  for  the  electrodes,  aluminum  for  the  current  collectors, 
and  a  salt  in  propylene  carbonate  or  acetonitrile  (or  an  equiva¬ 
lent  non-aqueous  solvent)  for  the  electrolyte  phase  [9,10,14]. 
A  large  effort  is  now  underway  in  the  electrochemical- 
materials  research  community  to  reduce  the  cost  of  the  acti¬ 
vated  carbon  [11-13],  the  remaining  high-cost  material  in 
these  systems.  Methods  by  which  one  can  assess  the  mi¬ 
crostructure  and  pore-size  distribution  of  the  activated  car¬ 
bons  relevant  to  capacitor  applications  are  reviewed  and 
discussed  in  [15].  In  addition,  because  of  the  high  spe¬ 
cific  power  and  power  density  of  these  capacitor  systems, 
it  may  be  possible  to  eliminate  the  dc/dc  converter  used 
in  propulsion  systems  for  controlling  voltages  between  en¬ 
ergy  storage  subsystems;  a  dc/dc  converter  is  used,  for  ex¬ 
ample,  when  a  battery  is  placed  in  parallel  with  a  fuel  cell 
[14]. 

This  paper  is  organized  as  follows.  Following  a  brief 
overview  of  the  experimental  tools  employed,  the  theoret¬ 
ical  underpinnings  and  governing  equations  utilized  in  the 
characterization  tests  and  adaptive  algorithm  are  presented. 
For  detailed  mathematical  modeling  of  electric-double-layer 
capacitor  systems  and  references  to  both  classical  and  re¬ 
cent  work  in  this  area  [16-25]  are  quite  helpful.  In  the 
last  portion  of  this  paper,  the  results  are  discussed  for 
the  characterization  studies,  the  use  of  the  current-voltage 
relations  as  an  adaptive  filter  is  examined,  and  finally 
the  adaptive  algorithm  is  implemented  to  post-process 
data. 


2.  Experimental 

The  capacitor  was  acquired  from  Saft  America.  It  has 
a  nominal  capacitance  of  3500  F  and  a  maximum  voltage 
of  2.8  V.  Electrochemical  tests  were  performed  on  an  Arbin 
BT2000  battery  testing  system  (Arbin  Instruments,  College 
Station,  TX).  Generally,  data  acquisition  rates  of  approxi¬ 
mately  1  point  per  100  ms  were  used.  Tests  were  conducted  at 
—30,  0,  25,  and  45  °C.  Temperature  control  was  provided  by 
a  Tenney  environment  chamber  (Model  Tenney20,  Lunaire, 
Williamsport,  PA).  At  each  temperature,  an  energy  test,  a 
charge  power  test  and  a  discharge  power  test  were  performed 
as  summarized  in  Table  1 .  Prior  to  a  test,  the  capacitor  was 
allowed  to  equilibrate  at  a  given  temperature  for  4  h. 


Table  1 

Characterization  test  procedures 

Energy  (SOC)  test 

1  Equilibrate  at  0  V 

2  Charge  at  100  A  to  1 .7  V 

3  Step  to  0  A  for  3  min;  record  Vo 

4  Discharge  at  1  A  until  0  V  is  obtained 

5  Repeat  step  2  in  increments  of  0. 1-2.4  V  (1.7,  1.8,  1.9,  . . .,  2.4  V) 

Discharge  power  test 

1  Equilibrate  at  0  V 

2  Charge  at  100  A  to  1 .7  V 

3  Step  voltage  to  1.4  V  for  5  s 

4  Step  to  0  A  for  3  min;  record  Vo 

5  Repeat  step  2  in  increments  of  0. 1-2.8  V  (1.7,  1.8,  1.9,  . . .,  2.4  V) 

Charge  power  test 

1  Equilibrate  at  0  V 

2  Charge  at  100  A  to  1 .7  V 

3  Step  voltage  to  2.8  V  for  5  s 

4  Step  to  0  A  for  3  min;  record  Vo 

5  Repeat  step  2  in  increments  of  0. 1-2.8  V  (1.7,  1.8,  1.9,  . . .,  2.4  V) 

These  capacitors  normally  operate  between  the  maximum  voltage  and  one- 
half  that  value.  The  maximum  voltage  for  this  system  is  2.8  V,  and  the  low¬ 
ermost  voltage  applied  was  1.4  V. 

2.7.  Energy  test 

The  capacitor  was  discharged  to  0  V  and  held  at  that  volt¬ 
age  for  300  s.  It  was  subsequently  charged  to  a  preset  potential 
of  1.7  V  at  a  constant  current  of  100  A  followed  by  an  open- 
circuit  rest  for  180  s.  The  capacitor  was  then  discharged  to 
0  V  at  a  constant  current  of  1  A.  This  sequence  was  repeated 
for  a  series  of  preset  potentials  at  0.1  V  increment  from  1.8 
to  2.4  V. 

2.2.  Discharge  power  test 

The  capacitor  was  discharged  to  0  V  and  held  at  that  volt¬ 
age  for  300  s.  It  was  then  charged  to  a  preset  potential  of 
1.7  V  at  a  constant  current  of  100  A.  Subsequently,  a  con¬ 
stant  voltage  of  1 .4  V  was  applied  for  5  s  and  the  capacitor 
was  allowed  to  rest  at  open-circuit  voltage  for  180  s.  This  se¬ 
quence  of  test  was  repeated  for  a  series  of  preset  potentials 
at  0. 1  V  increment  from  1 .8  to  2.4  V. 

2.3.  Charge  power  test 

This  test  was  similar  to  the  discharge  power  test,  except  a 
constant  voltage  of  2.8  V  was  applied  instead  of  1.4  V. 

3.  Equivalent  circuit  analysis  and  adaptive  algorithm 

The  model  used  for  the  characterization  analysis  and  the 
adaptive  algorithm  consists  of  two  parts.  There  is  an  electrical 
circuit  model  that  is  employed  to  describe  the  relationship 
between  the  currents  and  voltages  observed  at  the  terminals 
of  the  capacitor  (giving  rise  to  the  voltage-based  SOCy)  and  a 
coulomb-accumulation  model  that  describes  the  open-circuit 
voltage-based  on  the  history  of  currents  seen  by  the  capacitor 
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(giving  rise  to  the  coulomb-counting-based  SOCc),  which 
can  include  self  discharge  and  current  inefficiency  on  charge. 
The  electrical  circuit  model  is  illustrated  Fig.  1.  This  equiv¬ 
alent  circuit  is  used  to  extract  a  current-voltage  expression, 
which  is  then  employed  to  characterize  the  capacitor.  In 
addition,  the  current-voltage  expression  forms  the  basis  for 
the  voltage-based  SOCy.  It  should  be  recognized  that  even 
the  more  complex  equivalent  circuit  shown  in  the  upper 
schematic  of  Fig.  1  constitutes  a  greatly  simplified  approach 
to  modeling  a  capacitor  system  [25],  but  the  model  will 
prove  useful  for  the  purposes  of  an  adaptive  transversal 
filter. 

3.1.  Recursive  voltage  expression 

By  summing  the  voltage  around  the  circuit  shown  in  Fig.  1 , 
one  obtains: 


0 

V  =  —  +  IR, 
C 


(1) 


where  the  relationship  between  the  charge  on  the  capacitor 
and  the  associated  current  is: 


I  = 


d Q 

d  t 


The  initial  condition  for  V  or  I  corresponds  to: 

0(0) 

V(0)  =  ^  +  I(0)R. 

We  shall  assume  that  the  values  of  R  and  C  do  not  change 
substantially  over  the  time  scale  RC.  Differentiation  of  Eq. 
(1)  leads  to: 

dV  /  d  / 

—  =  -  +  R-. 
d  t  C  d  t 
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Fig.  1 .  Equivalent  circuit  representation.  The  upper  figure  provides  an  equivalent  circuit  for  a  capacitor  dominated  by  ohmic  drop  and  comprising  ideal  capacitive 
behavior  at  the  electrode-electrolyte  interface.  The  simplified  equivalent  circuit  used  to  construct  the  adaptive  filter  is  then  extracted  from  the  upper  figure  as 
indicated. 
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The  solution  to  this  equation  for  a  variable  current  source  is: 


0(0)  If 

V=^  +  IR+-  J  I |cd£ 

£=o 


For  constant  current  discharge, 


0(0) 

C 


+  IR  + 


It 

C' 


be  found  in  [26-32].  The  following  definitions  will  stream¬ 
line  notation: 

=  It  —  It— At  x2  =  (ft  + 

m\  =  R  m2  = 

y  =  ^meas,/  ^meas,/— A  t 

Eq.  (4)  can  now  be  recast  as: 
y  =  m\x\  +  m2v2, 


This  equation  is  used  to  extract  the  capacitance  versus  open- 
circuit  voltage  relationship  (Figs.  3  and  4).  Note  that  the 
open-circuit  voltage  Vo  corresponds  to  0/C  in  the  absence 
of  current  flow. 

The  first  step  in  constructing  the  voltage-based  portion  of 
the  SOC  algorithm  is  to  convert  the  integral  (2)  into  a  recur¬ 
sive  relation.  We  shall  work  through  the  first  few  time  steps  to 
elucidate  the  general  pattern  and  then  provide  the  necessary 
relation.  Initially,  t  =  t0  =  0,  and  Eq.  (2)  can  be  descretized  as: 


V\o  =  0lo^  +  /|otf. 

After  the  first  time  step,  t  =  t\,  and 

7|i  +  /|0  1 

V|i  =  (h  -  t0)  'U-  +  I\iR  +  (V|o  -  I\0R), 

where  the  integrand  I  has  been  approximated  by  the  average 
of  its  beginning  and  ending  values  of  the  time  step.  After  the 
second  time  step,  t  =  t2,  and 

V\2  =  fe  -  h)  -  +  I\2R  +  evil  -  I\iR), 


and  this  general  pattern  holds  for  time  step  N : 

I\n  +  7|^v_i  1 

v\N  =  (tN  -  tN- 1)  IA/  2  ^ 

+  I\nR  +  (V\N-i-I\N-iR). 


(4) 


Thus,  to  determine  the  voltage  at  time  A,  V\n,  one  only  needs 
the  present  value  of  the  current  and  the  previous  time-step 
values  for  the  current  and  voltage.  Consequently,  Eq.  (4)  is 
a  recursive  transformation  of  Eq.  (2)  into  a  descretized  time 
domain.  The  next  section  will  formulate  an  adaptive  proce¬ 
dure  for  the  estimation  of  R  and  C  from  a  history  of  currents 
and  voltages,  which  allows  for  the  state  estimation  of  the 
capacitor. 


3.2.  Weighted  recursive  least  square  for  the  estimation 
of  R  and  C 

Eq.  (4)  is  a  linear  relation  as  formulated,  and  least  squares 
methods  can  be  used  to  minimize  the  error  between  the  pre¬ 
dicted  voltage  V  and  the  measured  voltage  Vmeas  by  appropri¬ 
ate  adjustment  of  R  and  C  as  part  of  a  parameter  estimation 
scheme.  Reviews  on  the  substantially  interdependent  fields 
of  recursive  identification,  adaptive  filters,  optimal  estima¬ 
tors,  and  model-reference  adaptive  systems  (MRAS’s)  can 


and  the  error  s  that  we  seek  to  minimize  corresponds  to: 

N 

s  =  Wjlyj  -  (m\xij  +  m2X2,j)]2. 
j=  i 

Time  t  corresponds  to  index  j  =  N in  the  above.  Two  equations 
for  the  two  parameters  to  be  extracted  can  be  generated  by 
setting  the  partials  ds/dm\  and  ds! dm2  to  zero.  This  leads  to: 

1 

m  1  =  ^^(V\  yV2,2  ~  V2,yV\,2)  and 

=  ^-(V2,yVhl  -  V\  yV2  \),  (5) 

Det 

where  the  matrix  determinant  Det  is  given  by: 

Det  =  VL  IV2,2  ~  Vf. 

The  variances  corresponds  to: 


T/  ,  /  ,  1 

Vu,v\N  —  I  SU,v\ N  j  I  |  5 

\  Sw\n  /  ‘S’iiiIv 

in  which  the  sums  refer  to: 


N 

N 

sw  =  J2  Wj, 

su  =  £W 

j= 1 

7=1 

N 

N 

sv  =  J2  wivi . 

SU,v  =  ^2,  wjujvji 

j= 1 

7=1 

and  u  and  v  refer  to  x\,  jc2,  or  y.  It  should  be  noted  that 
the  matrix  system  is  symmetric,  Vp2  =  V2j.  Variances  have 
been  used  in  these  expressions,  along  with  the  normaliza¬ 
tion  associated  with  the  division  by  sw  in  order  to  keep 
the  resulting  matrix  elements  nearer  to  unity  (i.e.,  well 
scaled). 

Before  making  the  summations  recursive,  we  address  the 
weight  factor  wj.  There  are  two  reasons  why  one  would  want 
to  weight  data  sets  differently.  First,  some  observations  may 
be  subject  to  greater  disturbance;  here,  a  disturbance  refers 
to  a  phenomenon  not  accounted  for  in  the  system  model.  For 
example,  the  onset  of  secondary  reactions  during  capacitor 
charging  [16]  may  lead  one  to  discount  the  charge  obser¬ 
vations  relative  to  those  of  discharge,  as  the  impact  of  the 
secondary  reactions  is  not  treated  in  the  model  equations.  Sec¬ 
ond,  newer  observations  are  generally  more  important  than 
older  observations  in  determining  the  state  of  the  system, 
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Time,  s 

Fig.  2.  Energy  test  results  for  25  °C  data.  The  cell  is  charged  at  100  A  to  the 
indicated  maximum  voltage,  after  which  the  current  is  set  to  0  A  for  3  min, 
and  then  a  1  A  discharge  is  used  to  determine  the  coulombs  stored  in  the 
capacitor  at  the  equilibrium  voltage  Vo  (discharge  part  not  shown),  allowing 
for  the  determination  of  the  cell  capacitance  as  a  function  of  Vo  (cf.  Table  1 .) 
Analogous  plots  were  obtained  at  —30,  0,  and  45  °C. 

and  therefore,  should  usually  be  given  a  larger  weight  fac¬ 
tor  relative  to  older  observations.  For  these  two  reasons,  we 
decompose  the  weight  factor  into  a  time-weighting  factor  A 
and  a  general  weight  factor  y;  the  latter  of  which  can  be  used 


0  10  20  30  40  50  60  70  80 

Time,  s 


Fig.  4.  Cell  capacitance  determined  by  energy  test  (1  A  discharge)  as  a  func¬ 
tion  of  the  open-circuit  voltage  and  temperature. 

to  associate  a  specified  weighting  of  selected  events  (e.g., 
discharge  events  over  charge  events).  Hence, 

wj  =  ypN~j-  (?) 

It  can  be  shown  that  the  use  of  yields  an  exponential 
decay  in  the  influence  of  past  data  points  on  the  determination 
of  the  current  value  of  m\  and  m2 : 

XN~1  =  elnxN~j  =  e(N~j)lnX  «  f()r 

A  — >  1. 


Hence,  a  data  point  taken  100  s  prior  to  the  current  point 
has  less  than  40%  of  the  impact  on  the  regression  analysis 
relative  to  that  of  the  current  point.  A  data  point  taken  500  s 
prior  to  the  current  point  has  less  than  1  %  of  the  impact  on  the 
regression  relative  to  that  of  the  current  point.  More-complete 
discussions  on  exponential  forgetting  can  be  found  in  Section 
5.6.1  of  Ljnug  and  Soderstrom  [29],  Section  6.2  of  Anderson 
and  Moore  [31],  and  Sections  5.3  and  5.4  of  Kulhavy  [32]. 
The  summations  are  now  made  recursive  with  the  following 
definitions: 

Sw\n  =  Yn  +  Ksw\n-i), 

Su\n  =  YnUn  +  Ksu\n- l),  Sv\  N  =  YnVn  +  HsV\n~i), 

and 

Su,v\n  =  YNUNVN  +  k(sUjV\N-l), 
with  sUtV  =  sVtU.  Initially, 


Fig.  3.  Ohmic-corrected  voltage  versus  time  for  100  A  charge  (cf.  Fig.  2). 
The  variation  with  temperature  is  comparable  to  the  experimental  resolution. 
The  capacitance  fit  to  the  straight  line  corresponds  to  3410  F  (cf.  Eq.  (3)), 
which  is  consistent  with  the  1  A  discharge  data  shown  in  Fig.  4,  both  in 
terms  of  the  value  for  the  fit  capacitance  and  the  observed  invariance  of  the 
capacitance  with  temperature.  The  resistances  employed  to  form  V-IR  are 
plotted  in  Fig.  8. 


5yu  1 1  —  Y\  1 

su\i  =  Y\ui,  sv\i=yivi,  and.  (9) 

su,v\l  =  Yiuivi 

At  this  point,  the  equations  needed  to  regress  the  param¬ 
eters  m\  and  m2  have  been  fully  stated.  At  the  beginning  of 
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ated  with  the  energy  storage  device  (e.g.,  a  battery,  fuel  cell,  or 
capacitor)  to  ac  current  for  the  electric  machines  require  that 
Vmin  >0.5  Vmax  in  order  that  the  power  inverter  maintain  an 
acceptable  efficiency.  Thus,  for  the  figures  presented  herein, 
Vmin  is  set  to  zero  in  formulating  the  SOC,  and  we  investigate 
SOC’s  ranging  from  50  to  100  percent  in  the  power  capability 
analyses  (1.4  <  Vo  <  2.8  V);  100  percent  SOC  and  VmSLX  are 
taken  to  correspond  to  2.8  V,  which  implies  that  Vo  =  1.4  V 
refers  to  50%  SOC. 

We  can  also  calculate  an  SOC  based  on  coulomb  counting 
so  as  to  construct  a  current-based  SOC,  SOC/, 


SOC/(0  =  SO  C(t 


At)  + 


It  +  It- At 
(2  max  (2  min 


At 

T 


=  SOC(7 


1 

At)  -f  — 

c 


It  +  It- At 


Cm  ax  1 7=0  Vmin  |  /=() 


At 

T' 

(10) 


Fig.  5.  Cell  capacity  (Ah)  determined  by  energy  test  (1  A  discharge)  as  a 
function  of  the  open-circuit  voltage  and  temperature. 


each  time  step,  the  recursive  sums  (8)  are  calculated  (Eq.  (9) 
for  the  first  time  step),  which  are  then  substituted  in  the  vari¬ 
ances  (6);  upon  determining  all  variances,  the  parameters  m\ 
and  m2  are  deduced  by  means  of  Eq.  (5),  thereby  providing 
C  and  R. 


3.3.  State-of -charge  (SOC) 


The  previous  section  provided  for  the  extraction  of  C  and 
R  from  a  data  stream  provided  the  determinant  reflects  a 
well-posed  equation  system.  We  now  use  this  information 
along  with  coulomb  counting  to  construct  the  SOC  algo¬ 
rithm.  When  all  of  the  current  goes  to  charging  or  discharging 
the  electric-double-layer  associated  with  the  surface  of  the 
porous  electrodes  of  the  capacitor,  Eq.  (2)  can  be  rearranged 
such  that: 


2  =  2(0)  + 


s=t 

J  /|?d£  =  C(V 
/=0 


where  Q  is  the  magnitude  of  the  charge  on  the  (symmetric)  ca¬ 
pacitor  electrodes.  We  define  the  voltage-based  SOC,  SOCy, 
in  terms  of  the  minimum  charge  on  the  capacitor  2min  relative 
to  its  maximum  value  2max  allowed  under  normal  operating 
conditions, 

2  —  2 min  C(V  —  IR)  —  2min 

2  max  2min  CVmax|/=o  2min 

_  V  —  IR  —  Vmin  1 7=0 
Ejnax  I  / =0  Vmin  I  /=() 

where  2min  and  2max  correspond  to  the  minimum  and  max¬ 
imum  voltages  (Vmin  and  Vmax,  respectively)  under  zero  cur¬ 
rent  conditions.  In  vehicle  applications,  the  power  electronics 
(i.e.,  the  power  inverter)  used  to  convert  the  dc  current  associ- 


Both  the  voltage  and  current-based  SOC’s  contain  useful 
information,  and  a  weighted  average  is  thus  rendered  as  the 
final  (combined)  SOC, 

SOC  =  wSOCi  +  (1  -  w) SOCy,  (1 1) 

with  the  weight  factor  w  chosen  to  be  closer  to  1  for  enhanced 
stability  and  closer  to  zero  for  increased  responsiveness.  Un¬ 
less  otherwise  noted,  a  value  of  w  =  0.999  was  used  in  this 
work.  It  should  be  noted  that  w  must  be  viewed  in  the  context 
of  the  time  per  point,  as  the  coulomb  counting  includes  time 
integration,  unlike  the  more-responsive  voltage-based  SOC. 
Upon  inspection  of  the  formulas  for  SOCy  and  SOC/,  we 
note  that  SOCy  utilizes  the  adapted  resistance  R  and  SOC / 
relies  on  the  adapted  value  for  C,  and  in  this  context  both 
SOCy  and  SOC  j  are  adaptive  quantities  as  is  the  compos¬ 
ite  value  of  SOC.  In  addition,  it  is  important  to  note  that  in 
the  formulation  of  Eq.  (10)  for  SOC/,  we  increment  from 
the  previous  value  of  SOC  (not  SOC/),  thereby  linking  ex¬ 
plicitly  the  current-based  SOC  to  the  voltage-based  SOC. 
In  summary,  the  adaptive  algorithm  comprises  specifically 
Eqs.  (4)-(ll). 


3.4.  Power  capability 


The  maximum  discharge  power  can  be  expressed  as: 


Emax,  discharge  —  ^V  —  EVm in. 

That  is,  when  the  capacitor  voltage  obtains  its  lowest  ac¬ 
ceptable  value,  the  max  discharge  power  results.  First,  we 
consider  the  available  instantaneous  power;  i.e.,  the  power 
available  before  the  charge  on  the  capacitor  is  depleted  sig¬ 
nificantly  by  the  discharge  event.  In  this  case,  V  =  V(0)  +  IR, 
where  V(0)  is  the  system  voltage  at  zero  current  immediately 
prior  to  the  discharge  [V(0)  =  Q(0)/C,  cf.  Eq.  (2)],  and 


Em  ax,  discharge  —  IVmin  — 


[Vmin  ~  U(0)] 
R 


V 


mm' 


(12) 
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Fig.  6.  Interrupt  resistance  at  0  V  (upper  plot)  and  at  the  voltage  lid  (lower 
plot)  as  a  function  of  the  open-circuit  voltage  and  temperature. 


Similarly,  the  instantaneous  charge  power  for  the  ohmic  bat¬ 
tery  is  given  by: 


Anax,  charge  —  ^Emax  — 


[Vmax  ~  Vm 

R 


Vr 


max  ■ 


(13) 


The  ohmic  battery  does  not  address  transient  effects, 
which  are  important  for  times  that  are  greater  than  ~0.01 
RC.  In  this  case,  Eq.  (4)  can  be  inverted  to  solve  for  the  cur¬ 
rent,  and  consequently,  we  can  calculate  the  maximum  charge 
and  discharge  powers  available  for  the  time  interval  At: 


I\t  = 


Vt  -  [/,_Af  Ar/(2C)]  -  (V  -  IR)t_At 


R  +  [At/(2C)] 

Anax,  discharge  (  A  f)  =  ^Emin 

Emin  ~  [If- At  At/ (2C)]  ~  (V  ~  IR)t-At 
R  +  [At/(2C)] 
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Fig.  7.  Expanded  analysis  of  the  interrupt  resistance  at  the  voltage  lid.  The 
approximately  linear  trend  in  the  interrupt  resistance  versus  open-circuit 
voltage  depicted  in  Fig.  6  is  seen  to  apply  over  a  broad  voltage  range. 


and  Emax,  charge  (At)  —  IV max 

Vmax  -  Ut-At A//(2C)]  -  (V 


-  IR\- 


■At 


R  +  [At/(2C)\ 


V 


max  ■ 


(14) 


As  At  — >  0,  the  transient  power  calculations  match  those 
of  the  instantaneous  power  calculation.  To  implement  the 


Inverse  temperature,  1/(7",  K) 


Fig.  8.  Resistance  variation  with  absolute  temperature.  The  slope  yields  an 
activation  energy  of  2.0  kJ  mol-1.  The  resistance  values  correspond  to  an 
average  of  the  resistance  (at  the  indicated  temperature)  shown  in  the  upper 
plot  of  Fig.  6.  From  left  to  right,  the  temperatures  correspond  to  45,  25,  0, 
and  — 30°C. 
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Fig.  9.  Comparison  of  model  results  and  adaptive  filter.  Upper  plot:  model  results.  Lower  plot:  adaptive  filter.  Calculations  for  both  plots  employed  Eq.  (4). 
For  the  upper  plot,  V|yv_  i  in  Eq.  (4)  corresponds  to  the  previously  calculated  voltage,  whereas  V|yv_  i  corresponds  to  the  previously  measured  voltage  for  the 
lower  plot.  The  parameters  R  and  C  were  fit  by  least  squares  analysis,  with  each  point  receiving  equal  weighting. 


transient  equations,  the  respective  powers  are  calculated  im¬ 
mediately  after  the  algorithm  has  been  employed  to  finish 
the  SOC  determination  at  time  t.  In  this  case,  quantities  cal¬ 
culated  or  measured  at  time  t  are  then  stored  in  the  variables 
listed  in  the  respective  power  expressions  at  time  t  —  At.  Then 
one  must  state  the  duration  corresponding  to  the  desired  esti¬ 
mate  for  power.  For  example,  if  we  want  to  know  the  power 


estimates  2  s  from  “now”,  then  the  current  measured  and  ex¬ 
tracted  values  are  placed  in  the  t—  At  quantities,  t  and  At 
are  set  to  2  s,  and  the  right  sides  of  the  above  equations  yield 
the  desired  power  estimates.  Bounds  on  the  current  magni¬ 
tudes  should  also  be  observed  in  actual  applications;  that  is, 
the  power  capability  may  be  limited  by  a  maximum  current 
specified  for  the  capacitor  so  as  to  avoid  degradation. 
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Time,  s 

Fig.  10.  Expanded  view  of  plots  shown  in  Fig.  9.  The  nearly  linear  variation  in  the  cell  voltage  during  the  constant  current  charge  and  discharge  at  200  A 
magnitude  is  clearly  evident. 


4.  Results  and  discussion 

4.1.  Characterization  tests 

The  first  step  in  constructing  an  adaptive  algorithm  is  to 
characterize  the  system.  This  allows  one  to  determine  the 
suitability  of  selected  physical  models  (e.g.,  the  simplified 
equivalent  circuit  depicted  in  Fig.  1)  and  to  place  reasonable 


bounds  on  the  parameter  values.  This  latter  attribute  is  sig¬ 
nificant,  for  it  enables  one  to  determine  when  the  regressed 
parameters  take  on  reasonable  values,  indicating  robust  per¬ 
formance  of  the  algorithm  and  the  device.  The  three  types  of 
characterization  tests  implemented  in  the  work  are  described 
in  Table  1 . 

The  energy  test  can  be  used  to  measure  the  cell  capaci¬ 
tance  and  coulombic  capacity,  rendered  in  units  of  Ah  so  as 
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Fig.  11.  Charge  power  test  (cf.  Table  1).  The  extracted  state  of  charge  (SOC)  from  the  adaptive  algorithm  is  indicated,  as  is  the  voltage-based  contribution 
(SOCy)  and  that  corresponding  to  purely  coulomb  counting  (SOCc).  The  effect  of  the  specified  initial  SOC  and  the  algorithm  convergence  behavior  is  clarified 
by  comparing  the  upper  and  lower  plots;  an  erroneous  25%  initial  SOC  was  specified  in  the  lower  plot. 


to  be  consistent  with  battery  nomenclature.  Data  for  the  con¬ 
stant  current  (100  A)  charge  to  the  various  voltage  lids,  fol¬ 
lowed  by  the  180  s  open-circuit  stand,  are  provided  in  Fig.  2 
for  the  25  °C  condition.  As  indicated  in  Fig.  3,  the  effect 
of  temperature  on  the  data  as  plotted  was  difficult  to  dis¬ 
cern.  After  open-circuit  stand,  the  open-circuit  potential  Vo 
is  recorded.  It  then  remains  to  determine  how  much  charge 


resides  on  the  electrodes,  and  for  this  reason  a  low  current 
discharge  of  1  A  is  performed  to  0  V.  The  capacitance  C  is 
calculated  from  C  =  Q/V o,  where  Q  is  the  charge  passed  dur¬ 
ing  the  1  A  discharge.  The  measured  capacitance  C  and  ca¬ 
pacity  in  Ah  are  plotted  in  Figs.  4  and  5,  respectively,  over 
a  temperature  range  of  interest  for  automotive  applications. 
The  equilibrium  capacitance  plotted  in  Fig.  4  is  similar  to  the 
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Fig.  12.  Available  charge  power  analysis.  The  algorithm  utilizes  the  current  and  voltage  traces  shown  in  the  upper  plot  of  Fig.  1 1  to  produce  the  available  power 
results  depicted  above.  The  cell  power  is  the  current  multiplied  onto  the  voltage.  The  lower  plot  expands  the  abscissa  of  the  upper  plot  around  990  s.  The  key 
conclusion  is  that  the  available  power  predictions  rendered  by  the  recursive  algorithm  follow  closely  the  measured  power  through  the  current  reversal  at  984  s. 
For  this  analysis,  Z  =  0.995. 


effective  capacitance  extracted  from  the  slope  of  the  100  A 
charge  data  plotted  in  Fig.  3;  Eq.  (3),  with  Q  =  0  (correspond¬ 
ing  to  a  zero  initial  voltage  for  the  constant  current  charge), 
was  used  to  approximate  the  capacitance.  For  both  meth¬ 
ods  of  extracting  a  cell  capacitance,  little  temperature  vari¬ 
ation  is  observed.  The  curvature  in  the  voltage  versus  time 
curve  in  Fig.  3  is  consistent  with  porous  electrode  analyses 
[16]. 

When  the  current  is  changed  from  0  to  100  A  at  the  start 
of  the  experiment  corresponding  to  Fig.  2,  and  when  the  cur¬ 
rent  is  switched  from  100  A  to  zero  upon  reaching  the  volt¬ 
age  lid,  an  interrupt  resistance  can  be  ascertained.  Impedance 
due  to  ionic  diffusion  complicates  the  analysis;  for  our  pur¬ 
poses,  we  simply  plot  AV/AI  as  an  interrupt  resistance  in 
Fig.  6.  At  0  V,  the  interrupt  resistance  is  seen  to  be  nearly 
invariant  with  SOC  (i.e.,  Vo),  and  the  interrupt  resistance  de¬ 
creases  in  a  nearly  monotonic  manner  with  SOC,  a  feature 
more  fully  examined  in  Fig.  7.  A  succinct  explanation  for 


the  monotonic  trend  in  the  interrupt  resistance  seen  in  the 
lower  panel  of  Fig.  6,  and  in  Fig.  7,  remains  an  open  ques¬ 
tion.  When  the  logarithms  of  interrupt  resistances  near  0  V 
(values  averaged  over  the  SOC  for  each  temperature)  are  plot¬ 
ted  against  the  inverse  of  the  absolute  temperature  (Fig.  8), 
we  extract  a  rather  low  activation  energy  of  2  kJ  mol-1 .  This 
low  activation  energy  and  the  relative  invariance  in  the  ca¬ 
pacitance  with  temperature  help  to  explain  why  activated- 
carbon  electric-double-layer  capacitors  are  of  interest  for 
automotive  applications,  as  these  capacitors  maintain  their 
power  capability  at  lower  temperatures  much  more  so  than 
batteries. 

4.2.  Transversal  filter 

The  prior  characterization  of  the  capacitor  indicates  that 
the  RC  equivalent  circuit  depicted  in  Fig.  1  captures  many  of 
the  features  of  the  potential-current  relationship.  While  the 
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Fig.  13.  Extracted  (charge)  parameters  and  the  effect  of  the  forgetting  factor  A. 
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Fig.  14.  Discharge  power  test  (cf.  Table  1)  at  25  °C.  The  extracted  state  of  charge  (SOC)  from  the  adaptive  algorithm  is  indicated,  as  is  the  voltage-based 
contribution  (SOCy)  and  that  corresponding  to  purely  coulomb  counting  (SOCc)-  The  lower  plot  expands  the  time  axis  about  970  s  and  allows  one  to  view  the 
constant  current  (100  A  discharge)  prior  to  the  voltage  step  to  1.4  V,  the  minimum  voltage,  which  gives  rise  to  the  maximum  discharge  power  (IV). 


simplicity  of  the  RC  circuit  does  not  allow  one  to  understand 
the  underlying  phenomena  that  govern  the  potential-current 
relationship,  the  successful  correlation  of  the  capacitor  be¬ 
havior  does  bode  well  for  the  use  of  the  RC  circuit  for  the 
construction  of  an  adaptive  filter  and  ultimately  an  adaptive 
algorithm.  More  complex  models  that  include  the  numeri¬ 
cal  solution  to  coupled  partial  differential  equations  for  the 
purposes  of  describing  thermodynamic,  kinetic,  and  trans¬ 
port  phenomena  from  a  microscopic  basis  are  required  to 
improve  the  design  and  performance  of  activated-carbon 
electric-double-layer  capacitors.  A  comparison  of  the  RC  cir¬ 
cuit  as  a  model  versus  its  use  as  an  adaptive  filter  is  provided 


in  Fig.  9  (with  an  expanded  abscissa  shown  in  Fig.  10).  In 
both  plots,  the  measured  current  and  voltage  are  shown  by 
the  solid  lines.  The  current  source  was  constructed  so  as  to  be 
similar  that  which  would  be  used  in  a  capacitor-based  system 
in  an  automotive  application.  In  the  upper  plot  is  the  curve  de¬ 
picting  the  model  voltage  (symbols)  constructed  by  the  model 
equation  (cf.  Eq.  (4)): 


V\n  =  (Ov 


Ov-i) 


I\n  +  I\n  1 
2  C 


+  I\NR  +  (V\ 


calculated 
N- 1 


N- lR), 
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Fig.  15.  Available  discharge  power  analysis  for  25  °C  (upper  plot)  and  0°C  (lower  plot).  The  algorithm  utilizes  the  current  and  voltage  traces  shown  in  the 
upper  plot  of  Fig.  14  to  produce  the  available  discharge  power  results  depicted  above.  The  available  discharge  power  predictions  rendered  by  the  recursive 
algorithm  follow  closely  the  measured  power  through  the  voltage  reversals  (i.e.,  through  the  transition  from  a  constant  100  A  charge  to  a  1.4  V  step).  For  these 
analyses,  A  =  0.995.  Lower  discharge  powers  are  available  as  the  temperature  is  decreased,  consistent  with  the  results  of  Fig.  8. 


whereas  the  voltage  displayed  by  the  symbols  in  the  lower 
plot  corresponds  to  the  adaptive  filter, 


V\n  =  (?N 


tN- 1) 


I\n  +  ^Ia-i 
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C 


+  I\NR  +  (V\ 


measured 
N- 1 


I\n- lR). 


The  above  two  equations  differ  only  in  the  voltage  term 
VN~  i .  The  comparison  of  the  two  plots  makes  it  clear 

that  while  the  use  of  Eq.  (4)  is  not  appropriate  for  gen¬ 
eral  application,  the  implementation  of  Eq.  (4)  as  an  adap¬ 
tive  filter  does  provide  a  relatively  good  representation 
of  the  voltage  response,  and  the  extracted  parameters  de¬ 
duced  from  the  adaptive  filter  resemble  values  reported 


in  the  previous  section.  Thus,  as  is  commonly  the  case 
in  the  construction  of  model-reference  adaptive  systems, 
the  very  basic  model  provides  a  crude  approximation  to 
the  system  behavior,  while  the  recasting  of  the  model  as 
an  adaptive,  transversal  filter  [28]  provides  a  valuable  re¬ 
sult. 

4.3.  Adaptive  algorithm  results 

The  explicit  outputs  from  the  adaptive  algorithm  are  the  (i) 
state  of  charge,  (ii)  the  charge  and  discharge  power  capability, 
and  (iii)  the  regressed  parameters  (R  and  C,  which  can  be  used 
to  deduce  the  state  of  health  of  the  system).  Results  for  the 
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Fig.  16.  Extracted  (discharge)  parameters.  The  parameter  values  can  be  compared  with  those  depicted  in  Fig.  13  for  the  charge  analysis.  The  extracted  resistance 
R  is  larger  than  that  of  Fig.  13,  consistent  with  the  SOC  being  lower  for  the  discharge  tests  (Fig.  14,  about  50%  SOC  at  end  of  test)  than  for  the  charge  tests 
(Fig.  11,  about  90%  SOC  at  end  of  test);  the  same  trend  for  the  R  vs.  SOC  relationship  is  seen  in  the  lower  panel  of  Fig.  6.  Although  less  pronounced,  the  same 
explanation  applies  to  the  extracted  capacitance  C  (cf.,  Fig.  4). 


tests  labeled  “Discharge  power  test”  and  “Charge  power  test” 
in  Table  1  are  displayed  Figs.  11-16. 

The  performance  of  the  adaptive  algorithm  with  respect 
to  predicting  the  energy  left  in  the  system,  identified  with  the 
SOC,  is  depicted  in  Fig.  1 1.  In  the  upper  plot,  we  see  that  the 
algorithm  returns  a  value  close  to  that  obtained  using  sim¬ 
ply  coulomb  counting  (w  =  1  in  Eq.  (11)).  This  use  of  purely 
coulomb  counting  to  assess  the  accuracy  of  an  adaptive  al¬ 
gorithm  over  time  scales  long  relative  to  convergence  of  the 
algorithm  but  short  relative  to  the  influence  of  current  ineffi¬ 
ciencies  (side  reactions  [2,17])  has  been  shown  to  be  a  useful 
evaluation  method.  For  longer  durations,  as  in  actual  vehicle 
applications,  purely  coulomb  counting  will  become  substan¬ 
tially  inaccurate  due  to  current  inefficiencies.  The  duration 
for  convergence  for  the  chosen  set  of  parameters  is  indicated 
by  the  lower  plot  in  Fig.  1 1 ,  where  the  starting  SOC  for  the  al¬ 


gorithm  was  set  arbitrarily  to  25%,  instead  of  the  correct  value 
of  zero,  consistent  with  the  starting  state  of  equilibration  at 
0  V.  After  about  1400  s  the  voltage-based  SOC  is  nearly  iden¬ 
tical  to  that  of  the  combined  SOC.  For  this  analysis,  unless 
otherwise  noted,  w  =  0.999,  indicating  a  very  stable  algorithm 
greatly  influenced  by  SOC/.  The  algorithm  can  be  made  more 
adaptive  (and  less  stable)  by  lowering  w. 

The  prediction  of  the  capacitor’s  power  capability  asso¬ 
ciated  with  subjecting  the  discharging  (at  100  A)  capacitor 
to  an  abrupt  maximum  charge  condition  (2.8  V)  is  demon¬ 
strated  in  Fig.  12.  Generally  excellent  results  are  obtained, 
indicating  that  the  algorithm  and  the  power  capability  Eqs. 
(12)— (14)  work  quite  well.  The  lower  plot  in  Fig.  12  dis¬ 
plays  an  expanded  abscissa;  the  fact  that  the  solid  line,  re¬ 
flecting  the  actual  power,  tends  to  the  curve  represented  by 
the  symbols  (predicted  charge  power  capability)  when  the 
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potential  is  abruptly  set  to  2.8  V,  corresponding  to  a  current 
and  power  reversal,  implies  that  the  algorithm  correctly  cap¬ 
tures  the  maximum  charge  power  available  for  this  rather  se¬ 
vere  test  case.  The  charge  power  prediction  is  always  greater 
than  that  extrapolated  by  Eq.  (14)  when  At  =  500  ms  and 
is  close  to  the  instantaneous  charge  power  extrapolation, 
Eq.  (13). 

The  regressed  parameters  R  and  C  are  shown  in  Fig.  13. 
The  variation  in  the  extracted  parameters  can  be  reduced  by 
decreasing  the  forgetting  factor,  as  is  shown  in  the  lower  plot 
of  Fig.  13  for  the  0°C  condition;  however,  longer  conver¬ 
gence  times  are  then  expected.  The  general  lack  of  temper¬ 
ature  sensitivity  can  be  seen  by  comparing  the  0  and  25  °C 
results  in  Fig.  13;  because  the  results  tended  to  be  quite  sim¬ 
ilar  over  the  temperature  range  of  —30  to  45  °C,  in  terms 
of  algorithm  performance,  most  of  the  algorithm  analysis  is 
provided  for  the  25  °C  condition  only. 

The  analogous  experiments  to  investigate  the  discharge- 
power  prediction  capability  are  shown  in  Figs.  14,  15  and  16. 
Again,  0  and  25  °C  data  are  analyzed,  and  similar  values 
for  R  and  C  are  obtained.  As  with  the  test  of  the  charge- 
power  prediction  capability,  the  discharge  power  prediction 
is  always  greater  than  that  extrapolated  by  Eq.  (14)  when 
At  =  500  ms  and  is  close  to  the  instantaneous  discharge  power 
extrapolation,  Eq.  (12). 

5.  Summary  and  conclusions 

1 .  A  series  of  tests  ( cf.  Table  1 )  were  completed  to  charac¬ 
terize  electrochemically  the  performance  of  a  promising 
electric -double -layer  capacitor,  the  characterization  re¬ 
lies  on  a  simplified  equivalent  circuit  interpretation  ex¬ 
tracted  from  a  more-complete  mathematical  represen¬ 
tation  of  the  capacitor  system.  The  simplified  model 
amounts  to  a  series  RC  circuit,  and  the  resistance  R  and 
capacitance  C  are  extracted  over  the  temperature  range 
of  interest  for  automotive  applications  (—30  to  45  °C).  A 
variety  of  tests  yielded  consistent  values  for  R  and  C,  indi¬ 
cating  that  the  simple  model  is  sufficient  for  our  purposes. 
The  robustness  of  the  equivalent  circuit  enables  the  for¬ 
mulation  of  an  associated  model-reference  adaptive  algo¬ 
rithm:  the  measured  values  of  R  and  C  allow  one  to  place 
reasonable  bounds  on  the  adapted  values  of  these  param¬ 
eters  and  assess  whether  the  capacitor  system  is  function¬ 
ing  appropriately  in  a  vehicle  application  (i.e.,  deduce  the 
state  of  health  of  the  capacitor  system). 

2.  In  addition  to  the  high  power  capability  and  potentially 
low  costs  of  this  class  of  capacitors,  we  clarify  the  sub¬ 
stantial  invariance  of  the  device  performance  with  respect 
to  temperature,  a  distinct  advantage  over  battery  systems : 
however,  these  three  attributes,  relative  to  batteries,  come 
at  the  sacrifice  of  energy,  a  deficiency  that  renders  capac¬ 
itors  better  suited  to  supplement  fuel  cell  systems  rather 
than  internal  combustion  engines,  as  a  fuel  cell  can  be  the 
source  of  efficient  energy. 


3.  A  computationally  efficient  adaptive  algorithm  based  on 
weighted  recursive  least  squares  is  described  and  tested : 
the  algorithm  is  fully  recursive  in  that  the  only  variables 
required  for  on-line  regression  are  those  of  the  previous 
time  step  and  the  current  time  step.  A  time- weighting  tech¬ 
nique,  often  referred  to  as  exponential  forgetting,  is  em¬ 
ployed  to  damp  exponentially  the  influence  of  older  data 
points  on  the  regression  analysis.  Successful  comparisons 
of  the  algorithm’s  adapted  state  of  charge  and  power  pre¬ 
dictions  provide  an  initial  validation  of  the  algorithm. 
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